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Introduction 


This is a summary progress report of the NASA three- 
year Grant NAG 2-541 which was initiated in 
September 1988. This research project dealt with the 
development of new methodologies and algorithms 
for the multi-spectrum electromagnetic 
characterization of large scale non-metallic airborne 
vehicles and structures. As can be readily attested 
from the appended results of our research, a new 
robust, low memory and accurate methodology was 
developed which is particularly suited for modern 
machine architectures. This is a hybrid finite element 
method and combines two well known numerical 
solution approaches. That of the finite element 
method for modeling volumes and the boundary 
integral method which yields exact boundary 
conditions for terminating the finite element mesh. In 
addition, a variety of new high frequency results were 
generated (such as diffraction coefficients for 
impedance surfaces and material layers) and a new 
class of boundary conditions were developed which 
hold promise for more efficient simulations. During 
the course of this project nearly twenty five detailed 
research reports were generated(totalling more than 
1000 pages) along with an equal number of journal 
papers. The reports, papers and journals articles are 
listed in the appendices of this report along with their 
abstracts. The University of Michigan Radiation 
Laboratory report 02592 1-22-T is intended to provide 
a technical summary of the developed hybrid finite 
element method. A shorter version of this report has 
also appeared in the April 1991 issue of the IEEE 
Antennas and Propagation Transactions Magazine. 
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Below we summarize our progress on a year by year 
basis with recommendation on the continuation of this 
project. However, before doing so, a short statement 
of the project's initial goals and their importance to 
NASA and the Navy is given. 

Project Objectives 

• Implement a new formulation for scattering by two- 
dimensional(2D) and three-dimensional(3D) 
composite structures. At present, there is no available 
methodology for modeling large size non-metallic 
vehicles for electromagnetic characterization 
purposes. 

• New method is a modified version of the finite 
element method and claims a low O(n) memory 
demand making it attractive for large body 
simulations without any restrictions on the structure's 
composition. 

• Develop 2D and 3D general purpose codes and 
associated graphical user interfaces based on the new 
methodology. 

• Develop a new class of generalized impedance 
boundary conditions (GIBC) to allow the simulation 
of thin multilayer material coatings. These are to be 
incorporated with the finite element codes to model 
such coatings or to be used in place of absorbing 
boundary conditions(ABCs) 

• Develop diffraction coefficients for new geometries 
based on the proposed GIBCs. 

• Develop 2D and 3D numerical solutions based on the 
proposed GIBCs 


2 



Payoff 


• Proposed work will allow the simulation of three- 
dimensional composite vehicles and structures for test 
and evaluation. 

• To date, no other methodology can claim an O(n) 
memory demand without a comprise in accuracy and 
versatility. 

• No technique is presently available for treating 3D 
composite/penetrable structures 

• Because it provides near zone and far zone field 
distributions, it will also serve as a diagnostic tool in 
evaluating the effectiveness of surface treatments. 

• Complete, general purpose codes will be developed 
along with associated graphical user interfaces to be 
used in the test and evaluation of practical size 
vehicles. 


• The proposed generalized impedance boundary 
conditions should prove important in modeling high 
contrast RAM material either in the context of the 
proposed hybrid finite element method or in high 
frequency modelling. 
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Funding Profile 


YEAR AMOUNT 

PMTC. Univ.Mich. Total 

(Cost Sharing) 

1988- 89 $50,000 $29,000 $79,000 

1989- 90 $50,000 $29,000 $79,000 

1990- 91 $100,000 $25,000 $125,000 
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Project Milestones 


First Year Milestones 


• Develop 2-D formulation for the proposed hybrid 
finite element method (FE-CGFFT) 

• Implement and test methodology by comparison with 
results based on traditional approaches 

• Assess method's advantages and potential for future 
3D implementations 


• Develop higher order boundary conditions (GIBC) 
aimed at improving the simulation of multilayer 
coatings 

• Test accuracy and limitation of these boundary 
conditions 

• Begin development of new canonical analytical 
solutions by making use of the developed GIBC. 
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Second Year Milestones 


• Develop geometry generation package to be 
interfaced with the developed 2D hybrid finite 
element code. 

• Develop graphical user interface(GUI) for the 2D 
code, compatible with X-window platforms. 

• Start development of the proposed 3D hybrid finite 
element method (FE-CGFFT). 

• Apply FE-CGFFT method to a restricted class of 3D 
geometries; Test and evaluate method. 

• Develop analysis and begin code development for 
general Bodies of Revolutions (BOR) based on the 
new method. 


• Develop new diffraction coefficients based on the 
GIBC; 

Test and evaluate these for incorporation into GTD 
codes. 

• Develop 2D numerical solutions based on the derived 
GIBC. 
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Third Year Milestones 


• Complete and validate Body of Revolution Code and 
associated Graphical User Interface. 

• Develop edge-based shape functions and discretize 
equations using tetrahedral and other six-sided 
elements. 

• Develop test codes using the boundary integral or the 
ABCs for a simple structure such as a finite cylinder. 

• Develop and test other more accurate ABCs if 
necessary 

• Validate test codes by comparison with BOR moment 
method codes 

• Begin interfacing with 3D graphics packages such as 
SDRC IDEAS@ and PATRAN@. 


• Develop 3D numerical codes using simulations based 
on the Generalized Impedance Boundary 
Conditions(GIBCs) 

• Continue development of high frequency solutions 
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First Year Progress Summary 


• Developed analysis and associated computer code for 
2D scattering using the proposed finite element- 
boundary integral method in conjunction with the 
conjugate gradient-FFT (FE-CGFFT) 

• Wrote user interface and 2D mesh generation 
package. 

• Validated new code for rectangular structures and 
found its efficiency to be in accordance with original 
expectations. 


• Developed generalized impedance boundary 
conditions (GIBC) for multilayer material coatings on 
curved bodies. 

• Implemented and tested the developed GIBC for 
scattering by 2D coated bodies; A new numerical 
code was developed. 

• Generated a new canonical diffraction coefficient for 
skew incidence on impedance wedges. This 
diffraction coefficient will later prove to be of crucial 
importance in high frequency analysis of composite 
structures. 


• Five detailed reports were written on the above 
accomplishments. 

• Seven journal papers were submitted for publication 
and many conference publications resulted from this 
year's effort. 
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Second Year Progress 


• Developed formulation and code for scattering by 3D 
rectangular cavities, cracks and slots using the 
proposed FE-CGFFT method. Code demonstrated the 
method's unprecedented low memory, efficiency and 
versatility. 

• Made several improvements to the 2D FE-CGFFT 
code developed last year: 

• Developed new color graphics compatible with X- 
window platforms (post-processing tools). 

• Developed a new 2D mesh generator - AUTOMESH 
(pre-processing tool) 

•Validated code for a variety of composite geometries 
of arbitrary shape. 

•Determined that method's convergence can be 
improved by resorting to higher order boundary 
elements. 

• Developed the three dimensional FE-CGFFT analysis 
for Bodies of Revolution (BOR). 

• New BOR code was written and is being debugged. 

• Examined a number of commercial 3D mesh 
generation packages and found SDRC IDEAS® to be 
most useful for this method. 

• Developed new high frequency solutions for metal- 
dielectric junctions made possible by using the 
generalized impedance boundary conditions. 
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• Constructed a numerical solution based on the 
generalized boundary conditions for scattering by 2D 
cracks and gaps and impedance surfaces. 


• A total of fifteen(15) reports have been written during 
the first two years of this project. 

• More than twelve(12) journal papers have been 
published or submitted based on this work. 

• Many conference papers have also been presented. 


10 


Third Year Progress 


• Validated the Body of Revolution(BOR) Code and we 
are now improving its convergence characteristics by 
making use of higher order boundary elements. 

• Developed 3D finite element codes for non- 
rectangular and inhomogeneously filled cavities. 

• Developed new improved Absorbing Boundary 
Conditions (ABCs) to be used in conjunction with the 
proposed finite element solution. 

• Formulated new edge-based tetrahedral elements 
suitable for electromagnetic simulations. Coded and 
tested formulation for modeling closed cavities 

• The development of a general purpose 3D finite 
element - ABC code using tetrahedral elements is 
nearly complete. 


• New diffraction coefficients were generated for 
dielectric-to-dielectric junctions using the developed 
GIBCs. 

• The non-uniqueness of GIBCs was investigated and 
resolved for a number of cases. 

• Developed solutions for creeping and surface waves 
on multilayer dielectrics. A code was also written and 
tested. 

• Developed a 3D formulation and code for multilayer 
dielectric coatings using the developed GIBCs. 
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• Twenty two (22) reports totalling more than 1200 
pages have been written so far as a result of this 
effort. 

• Over thirty journal papers have been written to date. 

• Over twenty conference papers have also been written 
and presented. 
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Summary of Accomplishments 


During the course of this study we have developed 
and validated a new method for scattering by 
composite vehicles. The method has proven superior 
to existing approaches by requiring a low O(n) 
storage without compromising the accuracy of the 
solution. So far, up to 150,000 unknown systems 
have been solved via the method with rather 
remarkable convergence rates. 

Based on this method, a general purpose scattering 
code was written along with a mesh generator 
compatible with X window platforms. A Body of 
Revolution code was also written as a first step in 
moving toward a full scale three dimensional 
implementation. At this point the three-dimensional 
finite element formulation has been developed and 
coded for edge-based elements which are best suited 
for electromagnetic simulations. A three dimensional 
code based on a new class of absorbing boundary 
conditions should be ready within one to two months. 
The accuracy of this code will be examined before 
proceeding with the incorporation of the exact 
boundary integral equations for terminating the mesh. 


In addition to the development of the new numerical 
method, a rather successful study was carried out 
pertaining to the development of improved 
impedance boundary conditions and diffraction 
coefficients for impedance and material edges. These 
were crucial in improving the capability of existing 
high frequency codes. In addition, they provide an 
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efficient and elegant simulation of thin layers and 
impedance surfaces in the context of finite element 
solutions. 

Future Work 


During the previous grant we developed and tested 
the proposed hybrid finite element method (FE- 
CGFFT). The method was applied to a variety of 2D 
and some 3D applications but by and large this initial 
study was concerned with the development and 
validation of the method. However, the full potential 
of the method remains to be explored and this should 
be the driving goal of the next effort. It is essential 
that the full 3D implementation of the method be 
completed as planned and be applied for a 
characterization of new large scale composite 
geometries. New discretization elements should be 
examined and implemented in a manner that leads to 
improved system conditions. Further, preconditioning 
schemes should be examined and new solvers for 
vector and multi-processor architectures be explored 
for solving very large systems involving millions of 
unknowns. A major advantage of the proposed 
methodology is its suitability for solution on multi- 
processor architectures and this feature should be 
explored to its maximum. Finally, the role of pre- and 
post-processing software should not be ignored. 
When dealing with large structures, the generation of 
geometry discretization files is a major, time 
consuming task and must be well integrated with the 
analysis code. Also, the processing and presentation 
of the multi-spectral RCS or near field data present us 
with a challenging task and ways must be developed 
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to effectively display the information contained by 
such data. 
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PRESENTATION VIEWGRAPHS 
ON THE 

FORMULATION, VALIDATION AND 
APPLICATION OF THE 
HYBRID FINITE ELEMENT-CGFFT METHOD 
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Geometries where the proposed Finite 
Element-Boundary Element has been Applied 
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May converge slowly unless 
conditioned. 



Finite Element-Boundary Integral 
CGFFT Hybrid Method 
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Formulation of the Proposed Finite 
Element-Boundary Element Method 


z 




• Tightly enclose structure by a contour (for 2D) or a surface 
(for 3D). This enclosure separates the infinite space into - 
region Q and Qoo. 

• Formulate and Discretize the field in the interior region Q via 
the finite element method. 

- use triangles for 2D 

- use tetrahedra or bricks for 3D 

• Employ the Boundary integral equation on T to relate the 
tangential fields. Discretize this i.e. via the Boundary 
Element method. 

• Combine the Finite element and boundary element matricies 
to solve for the modal or element fields. 

• If r is kept flat or circular, the boundary integrals are 
convolutional. 

• By using an iterative solution algorithm to solve the system, 
the FFT can be employed to evaluate the convolutional 
boundary integrals leading to an O(n) storage requirement. 
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Finite Element Formulation 
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Linear interpolation: 


Assume 

4> e ( x ,y) 

Then 



3 

<!> e (x,y) = ^ Nj(x,y)4>l 

i=l 

where 


Nf(x, y ) = + b { x + Ciy) 

and 

a l = X 2V3 - V2 X 3 I &1 = 2/2 “ 2/3 i c l = x 3 ~ x 2 

a 2 ~ x 3Vl ~ V3 X 1 I &2 = 2/3 — 2/1 i c 2 = X 1 ~ x 3 

a 3 — x l2/2 ~ Vl x 2 I ^3 = 2/1 — 2/2 I c 3 = x 2 — ^1 



1 x 1 yi 
1 x 2 2/2 
1 *3 2/3 


= area of the eth element 
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Discretization or tne Functional 
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jk 0 Z 0 f L*L*j I <n 


Boundary Integral Formulation 
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Backscatter Echo Width (dB) 




0.0 30.0 < 50.0 90.0 120.0 150.0 180.0 

Angle (deg) 


Figure 1 : E z and H z backscatter echowidth patterns for the illustrated 

perfectly conducting cylinder. 
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FILE: missile_out 
CREATION DATE: 1 Aug 1990 
STRUCTURE: missile shape 
ENCLOSURE: rectangle 

number of nodes: 636 

number of elements: 90S 

nodes on pec boundary: 1 72 

nodes on obs boundary: 1 92 

nodes on unknowns (E-pol): 464 
nodes on unknowns (H-pol): 636 


QQ 

T3 



flQ 


Missile Structure 



Figure 7: Backscatter patterns for a missile-like perfectly conducting and 

coated cylinder. 
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COATED OGIVE 


p = S . it 



Figure 2(b): Example of a rectangular mesh enclosing a coated ogive. 


FILE: ogive_Ieo_out 
CREATION DATE: 7 Aug 1990 
STRUCTURE: 0.5 x 1 X coated conductor 
ENCLOSURE: rectangle 


number of nodes: 1084 

number of elements: 1936 

nodes on pec boundary: 140 

nodes on obs boundary: 92 


nodes on unknowns (E-poi): 944 
nodes on unknowns (H-pol): 1084 



15.0 30.0 45.0 60.0 


Angle [degj 


Figure 2(c): H-polarization backscatter pattern for the coated ogive in 

Fig. 2(b). 
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Coated Trailing Edge FEM Mesh 


case 2 


±W l-(x/2.5) 2 

y. = 2 

m (±0.8232 A(x) 

-2.5<x<0 

0<x<2.5 

A(x) = V(l-(x/2.7182pj 

- 0.3926 

y out = +0.8116 B(x) 

0.1<x<3.0 


B(x) = V l-[(x-.l)/3.1416] 2 -.3846 


£ r =^-j l 



5 } 
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Backscatter Echo Width [dB/X] Backscatter Echo Width [dB/X] 


Case 2: Coated trailing edge, (£=2-jl), H-pol 



Angle [deg] 


Case 2: Coated trailing edge, (£=2-jl), E-pol 
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with a material having c r — 2.56 — j0.256, — 1-2 jO.l, tn 1A. d 



Top and bottom layers: e r = 2.56, n T = 1-0, 0.2 cm thick ^.middle layer: e T - 4 .0. 


H r = 1.0, 0.2 cm thick. 
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2cm 
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('on'* v/d 



c = 0.25 A. material filling 

a=lX e r = 7 - j0.5 ( a ) = <)> = 90° 

b = 0.25A, jx r = 1-8 - j.01 . (b) = 0 = 0° 



0 . 10 . 20 . 30 . 40 . 50 . 60 . 70 . 80 . 90 . 0 . 10 . 20 . 30 . 40 . 50 . 60 . 70 . 80 . 90 . 

0 (degrees) 0 (degrees) 

(a) (b) 


Backscatter patterns for a material-filled cavity versus incidence angle, a 
1.0A, b = 0.25A, c = 0.25A, c r = 7.0 - j0.5, Mr = 1.8 - jO.l. Solid and dashed lines 
represent the solution of this method; circles and squares represent the moment 
method/modal solution, (a) 9 = 90°. (b) 9 = 0°- 
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10 . 

0. 

10. 

20. 

30. 

40. 

50. 

60. 



: 0 (degrees) 


0 (degrees) 

(b) 


Backscatter patterns for a multilayer material-filled cavity versus incidence an- 
gle; a = 1.0A, b = 0.25A, c = 0.25A. Top layer: e T = 7.0 — j’0.5, = 1.8 — j'0.1. 

0.0625A thick; middle layer: c r = 3.0 — j'0.05, /x r = 2.0, 0.125A thick; bottom 
layer: c r = 5.0 — jO.03, n r — 1.0, 0.0625A thick. Solid and dashed lines represent 
the three-dimensional solution; circles and squares represent the two-dimensional 
solution, (a) </> = 90°. (b) if = 0°. 
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Finite Element Formulation 

— 

t. 

z 

• Subdivide the region internal to r a into finite elements 

• Expand the fields over element e as 

= R £ eUi X;(R, Z) 
i=i 

= R X Nj(R, Z) 

j = 1 

• Apply the method of weighted residuals to the CAP equa- 
tions over the eth element 

JJ S , R N‘(R,Z ) R e sl dS = 0 
JJ S e R N*(R,Z) R.' s2 dS = 0 

where 

R e sl = V • [/ m (e r i?V^ + x V^)] + ^ 

^2 = V • [/ m (/i, x Vtf)] + ^ 

and A^- e is the usual linear shape function. 
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FE-BE System Assembly 


• The final finite element - boundary element system is 


A aa A a I B aa 
Ala An 0 
Lnn 0 M, 


aa 



^ m<t>a 


0 


^ m<f>I 

— 

0 


'U'mt 


me 

u m<pa 


where 


Paa 


P -Q 
Q P 


M-aa 


-P l -Q l 
Q l —P l 


and 


p = m - p*]c 
Q = IQ*C + Q*'D 


A typical matrix element is 


Qt 2 ].. = (R2,R2,Z i+i ,Z')Z'dZ' 

• Solve the equation via the CG algorithm and use the FFT 
to compute some of the integrals 
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Conducting Sphere (p=0.5A,) 









lap] 
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Frequency (GHz) 
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Finite Element Boundary Integral Method 


METHOD' S EFFICTFNrV; 

( 1 ) Expected to be very efficient in 
computation time (we have observed 
unusually fast convergences) 

(2) Less efficient for RCS vs angle 
computations 

( 3 ) Very efficient in memory use: 

Memory demand is approx. 

5 Nt + 10 Ns complex 
5 Nt + 2 Ns integers 
Where 

Nt = total number of unknowgtfand 
Ns = number of aperture unknowns 
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Finite Element Boundary Integral Method 


METHOD'S VERSATILITY AND 
CAPABILITY: 

(1) Treats layered substrates and 
superstrates 

(2) Treats microstrip patch antennas and 
slot antennas 

(3) Models impedance loads and short- 
circuit pins 

(4) Models probe feeds, microstrip line 
edge feeds and microstrip sandwich 
feeds 

(5) Provides characterization of both 
scattering and radiation properties 

(6) Has O(n) memory demand and is 
capable of treating large geometries 
with tens of thousands unknowns 
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Report no. 025921- 1-T 


Abstracts of technical reports 


DERIVATION AND APPLICATION OF A CLASS OF GENERALIZED 
BOUNDARY CONDITIONS - II 

J. L. Volakis, T. B. A. Senior, and J.-M. Jin 
The Radiation Laboratory 

Department of Electrical Engineering and Computer Science 
The University of Michigan 
Ann Arbor, Ml 48109-2122 


Abstract 


Boundary conditions involving higher order derivatives are presented for 
simulating surfaces whose reflection coefficients are known analytically, numerically or 
experimentally. Procedures for determining the coefficients of the derivatives are 
discussed, along with the effect of displacing the surface where the boundary 
conditions are applied. Provided the coefficients satisfy a duality relation, equivalent 
forms of the boundary conditions involving tangential field components are deduced, 
and these provide the natural extension to non-planar surfaces. As an illustration, the 
simulation of metal-backed uniform and three layer dielectric coatings is given. It is 
shown that fourth order conditions are capable of providing an accurate simulation for 
the uniform coating at least a quarter of a wavelength in thickness. Provided, though, 
some comprise in accuracy is acceptable, it is also shown that a third order condition 
may be sufficient for practical purposes when simulating uniform coatings. 
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Scattering by a Two Dimensional Groove in a Ground Plane 

Kasra Barkeshli and John L. Volakis 


Abstract 

Higher order boundary conditions involve derivatives of the fields beyond the first 
and were recently shown to be more effective than traditional first order conditions in 
modeling dielectric coatings and layers. In this report an application of a third order 
generalized boundary condition to scattering by a filled rectangular groove is 
presented. Deficiencies of such higher order boundary conditions are addressed and 
a correction is proposed for the present case. As part of the process of examining and 
Improving the accuracy of the proposed generalized boundary conditions, an exact 
solution is developed and a comparison is provided with a solution based on the 
standard impedance boundary condition. 
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Report no. 02592 1-3-T 


Preface 

This semi-annual report describes our progress during the period 
from September 1988 to February 1989. There are four separate 
tasks described in this report. Although, some of these tasks are 
interdependent, their reports should be read independently. That 
is, figure and reference numbering is consecutive only within the 
description of a given task. As can be expected, the progress report 
on each task is brief. The reader is expected to refer to the detailed 
reports or papers referenced in the progress reports. 
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Report no. 02592 1-4-T 


An Approximate Skew Incidence Diffraction Coefficient for an 

Impedance Wedge 


Hasnain H. Syed and John L. Voiakis 
Radiation Laboratory 

Department of Electrical Engineering and Computer Science 
The University of Michigan 
Ann Arbor, MI 48109-2122 


Abstract 

An approximate diffraction coefficient is presented for an impedance wedge at skew 
incidence. This is based on suitable modifications of the available skew incidence solution 
for an impedance half plane. The derived diffraction coefficient reduces to the known ones 
for an impedance wedge at normal incidence, a perfectly conducting wedge at skew 
incidence as well as an impedance half plane at skew incidence. An additional crucial 
requirement that the total field remains continuous at die reflection and shadow boundary is 
also imposed. Using the derived approximate diffraction coefficient a first order solution is 
obtained for the scattering by a polygonal impedance cylinder at skew incidence. This is 
then compared with numerical data which validate the accuracy of the presented dyadic 
coefficient. 
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Report no. 02592 1-5-T 


perivation of Generalized Transit ion /Boundary Conditions 
for Planar Multiple Layer Structures 1 

September, 1989 

M. A. Ricoy and J. L. Volakis 
Radiation Laboratory 

Department of Electrical Engineering and Computer Science 
The University of Michigan 
Ann Arbor, Michigan 48109-2122 

i4fcs<rac<-Infinite order generalized impedance boundary conditions (GIBCs) and 
generalized sheet transition conditions (GSTCs) for planar multilayer configura- 
tions are developed via the Taylor series expansion method. The conditions are 
derived in a matrix product form where each matrix corresponds to a specific 
layer. An overall composite boundary/ transition condition is obtained by mak- 
ing finite order approximations to the elements of each matrix for the cases of 
“low contrast” and “high contrast” material layers. The accuracy of the trun- 
cated boundary conditions is examined by comparing their implied reflection and 
transmission coefficients with the corresponding exact coefficients. Design curves 
are also given which relate the maximum order of the conditions required to sim- 
ulate a coating or layer of specific thickness and contrast. Expressions are then 
derived for the reflection and transmission coefficients of the GIBC/GSTC sheets 
and these are compared to exact coefficients to demonstrate the validity of the 
derived GlBCs/GSTCs. 

'This work *u supported in part by the NASA Ames Research Center under Grant NAG 2-MI. 
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Report no. 02592 1-6-T 


A Combined Finite Element-Boundary Element Formulation 
for Solution via CGFFT of 2-Dimensional Scattering Problems 

J.D. Collins and J.L. Volakis 


Radiation Laboratory 

Department of Electrical Engineering and Computer Science 
The University of Michigan 
Ann Arbor, MI 48109-2122 


Abstract 

A new technique is presented for computing the scattering by two-dimensional structures 
of arbitrary composition. The proposed solution approach combines the usual finite ele- 
ment method with the boundary integral equation to formulate a discrete system. This 
subsequently solved via the conjugate gradient (CG) algorithm. A particular characteris- 
tic of the method is the use of rectangular boundaries to enclose the scatterer. Several of 
the resulting boundary integrals are therefore convolutions and may be evaluated via the 
fast Fourier transform (FFT)in the implementation of the CG algorithm. The solution 
approach offers the principle advantage of having 0(N ) memory demand and employs 
& one- dim ensional FFT versus a two-dimensional FFT as required with a traditional 
implementation of the CGFFT algorithm The speed of the proposed solution method 
is compared with that of the traditional CGFFT algorithm, and results for rectangular 
bodies are given and shown to be in excellent agreement with the moment method. 
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Preface 


This semi-annual report describes our progress during the period from February 
1989 to September 1989. In all, six technical reports, seven journal articles and eight 
conference papers have been completed since the beginning of this grant one year ago. 
These are referenced at the end of each task described in the reprint. 

There are four tasks described in this report. Each should be read 
independently. That is, figure and reference numbering is consecutive only within the 
description of the task. As can be expected, the progress reports are very brief and the 
reader should to refer to the referenced technical reports for detailed coverage. 
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Report no. 02592 1-8-T 


Diffraction by a Multilayer Slab Recessed in a Ground Plane 
via Generalized Impedance Boundary Conditions 1 

December, 1989 

M. A. Rjcoy and J. L. Volakis 
Radiation Laboratory 

Department of Electrical Engineering and Computer Science 
The University of Michigan 
Ann Arbor, Michigan 48109-2122 

>46s/ract-The diffraction problem associated with a multilayer material slab re- 
cessed in a perfectly conducting ground plane is formulated and solved via the 
Generalized Scattering Matrix Formulation (GSMF) in conjunction with the dual 
integral equation approach. The multilayer slab is replaced by a surface obeying a 
generalized impedance boundary condition (GIBC) to facilitate the computation 
of the pertinent Wiener Hopf split functions and their zeros. Both E z and H z 
polarizations are considered and a number of scattering patterns are presented, 

some of which are compared to exact results available for a homogeneous recessed 
slab. 
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Report no. 02592 1-9-T 


An Asymptotic Analysis of the Plane Wave Scattering by a 
Smooth Convex Impedance Cylinder 


Hasnain H. Syed and John L. Volakis 
Radiation Laboratory 

Department of Electrical Engineering and Computer Science 
The University of Michigan 
Ann Arbor, MI 48109-2122 


Abstract 

A rigorous UTD (Uniform Geometrical Theory of Diffraction) analysis of the 
diffraction by a smooth convex impedance cylinder is presented. The analysis parallels the 
one employed for the perfectly conducting cylinder. Ray solutions are obtained which 
remain valid in the transition region and uniformly reduce to those in the lit and the shadow 
regions. In addition, a ray solution is presented when the observation point is in the close 
vicinity of the cylinder. The resulting transition region expressions are in terms of Fock- 
type integrals which are evaluated using a numerical scheme based on Fourier quadrature. 
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A prNlTE ELEMENT-BOUNDARY INTEGRAL FORMULATION 
FOR SCATTERING BY THREE-DIMENSIONAL 
CAVITY-BACKED APERTURES 

Jian-Ming Jin and John L. Volakis 
Radiation Laboratory 

Department of Electrical Engineering and Computer Science 
The University of Michigan 
Ann Arbor , Michigan 48109 

ABSTRACT 

A new numerical technique is proposed for the electromagnetic characterization of 
the scattering by a three-dimensional cavity-backed aperture in an infinite ground plane. 
The technique combines the finite element and boundary integral methods to formulate 
a system of equations for the solution of the aperture fields and those inside the cavity. 
Specifically, the finite element method is employed to formulate the fields in the cavity 
region and the boundary integral approach is used in conjunction with the equivalence 
principle to represent the fields above the ground plane. Unlike traditional approaches, 
the proposed technique does not require knowledge of the cavity’s Green’s function and is, 
therefore, applicable to arbitrary shape depressions and material fillings. Furthermore, 
the proposed formulation leads to a system having a partly full and partly sparse as well 
** symmetric and banded matrix which can be solved efficiently using special algorithms. 
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Abstract 

A new method for the computation of electromagnetic scattering from arbitrary two- 
dimensional bodies is presented. The method combines the finite element and boundary 
clement methods leading to a system for solution via the Conjugate Gradient FFT al- 
gorithm. Two forms of boundaries aimed at reducing the storage requirement of the 
boundary integral are investigated. It is shown that the boundary integral becomes con- 
volutional when a circular enclosure is chosen, resulting in reduced storage requirement 
when the system is solved via the Conjugate Gradient FFT method. The same holds for 
the ogival enclosure, except that some of the boundary integrals are not convolutional 
and must be carefully treated to maintain the O(N) memory requirement. Results for 
several circular and ogival structures are presented and shown to be in excellent agree- 
ment with those obtained by traditional methods. 
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Preface 


io»Q m our progress during the period from September 

1989 to February 1990. Several technical reports and papers have been written and these 
are listed at the end of each task. 

There are four tasks described in this report. Each should be read independently 
That is figure and reference numbering is consecutive only within the description of the 
task. As can be expected, the progress reports are very brief and the reader should refer to 
the referenced technical reports for detailed coverage. 
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Preface 


This semi-annual report describes our progress during the period from February 
1990 to September 1990. Several technical reports and papers have been written and these 
are listed at the end of each task. 

There are two tasks described in this report. Each should be read independently. 
That is, figure and reference numbering is consecutive only within the description of the 
task. As can be expected, the progress reports are very brief and the reader should refer to 
the referenced technical reports for detailed coverage. A total of sixteen technical reports 
have already been submitted and two more are currendy being prepared. Also, more than 
18 publications to refereed journals and more than 17 conference papers have resulted from 
this sponsored research. 
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Electromagnetic Scattering From Two-Dimensional Thick 

Material Junctions 


August, 1990 


M. A. Ricoy and J. L. Volakis 
Radiation Laboratory 

Department of Electrical Engineering and Computer Science 
The University of Michigan 
Ann Arbor, Michigan 48109-2122 


/4&5*rac<- Because material junctions are commonplace on structures whose radar cross section is 
of interest, it is essential that their scattering properties be adequately characterized The stan- 
dard impedance boundary condition (SIBC) has been employed in the past along with function 
theoretic techniques to develop simple scattering models of material junctions with thin and/or 
high loss slabs. To extend these models to more general slabs, generalized impedance boundary 
conditions (GIBCs) and generalized sheet transition conditions (GSTCs) have been proposed. 
Unfortunately, the solutions obtained with these are usually non-unique in the form of unknown 
constants, and although the constants have been resolved for a few special cases, previous efforts 
were unable to determine them in the general case. 

This report examines the problem of the plane wave diffraction by an arbitrary symmetric 
two-dimensional junction , where Generalized Impedance Boundary Conditions (GIBCs) and Gen- 
eralized Sheet Transition Conditions (GSTCs) are employed to simulate the slabs. In chapter 2, 
GIBCs and GSTCs are constructed for multilayer planar slabs of arbitrary thickness and the 
resulting GIBC/GSTC reflection coefficients are compared with exact counterparts to evaluate 
the GIBCs/GSTCs. In chapter 3 the plane wave diffraction by a multilayer material slab recessed 
in a perfectly conducting ground plane is formulated and solved via the Generalized Scattering 
Matrix Formulation (GSMF) in conjunction with the dual integral equation approach. Various 
scattering patterns are computed and validated with exact results where possible. 

In chapter 4, the diffraction by a material discontinuity in a thick dielect ric/ferrite slab is 
considered by modelling the constituent slabs with GSTCs. A non-unique solution in terms of 
unknown constants is obtained, and these constants are evaluated for the recessed slab geometry 
of chapter 3 by comparison with the solution obtained therein. Several other simplified cases 
w also presented and discussed. In chapter 5 an eigenfunction expansion method is introduced 
to determine the unknown solution constants in the general case. This procedure is applied to 
the solution of chapter 4, and scattering patterns are presented for various slab junctions and 
compared with alternative results where possible. Chapter six presents a short summary of this 
report and some recommendations for future work. 
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High Frequency Scattering by a Smooth 
Coated Cylinder Simulated with 
Generalized Impedance Boundary 
Conditions 

Hasnain H. Syed and John L. Volakis 


Radiation Laboratory 

Department of Electrical Engineering and Computer Science 
The University of Michigan 
Ann Arbor, Michigan 48109-2122 


Abstract 

Rigorous UTD (Uniform Geometrical Theory of Diffraction) diffrac- 
tion coefficients are presented for a coated convex cylinder simulated 
with generalized impedance boundary conditions. In particular, ray 
solutions are obtained which remain valid in the transition region and 
reduce uniformly to those in the deep lit and shadow regions. These 
involve new transition functions in place of the usual Fock-type inte- 
grals, characteristic to the impedance cylinder. A uniform asymptotic 
solution is also presented for observations in the close vicinity of the 
cylinder. As usual, the diffraction coefficients for the convex cylinder 
are obtained via a generalization of the correpsonding ones for the 
circular cylinder. 
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Abstract 

A new method is presented for the computation of electromagnetic scattering from 
axially symmetric bodies. To allow the simulation of inhomogeneous cross sections, the 
method combines the finite element and boundary element techniques. Interior to a 
fictitious surface enclosing the scattering body, the finite element method is used which 
results in a sparse submatrix, whereas along the enclosure the Stratton- Chu integral 
equation is enforced. By choosing the fictitious enclosure to be a right circular cylinder, 
most of the resulting boundary integrals are convolutional and may therefore be evaluated 
via the FFT with the system is iteratively solved. In view of the sparse matrix associated 
with the interior fields, this reduces the storage requirement of the entire system to O(N) 
making the method attractive for large scale computations. This report describes the 
details of the corresponding formulation and its numerical implementation. 
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Preface 

This semi-annual report describes our progress during the 
period from September 1990 to February 1991. Two specific tasks 
are described and each should be read independently. That is, figure 
and reference numbering is consecutive only within the description 
of the task. As can be expected the progress reports are very brief 
and the reader should refer to the referenced technical reports for 
detailed coverage. A total of 23 journal articles, 17 conference 
publications and 18 technical reports have been written since the 
beginning of the Grant. A complete list of these is given in pp. 40- 
45. 
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A FINITE ELEMENT - BOUNDARY INTEGRAL METHOD FOR 
SCATTERING AND RADIATION BY TWO- AND 
THREE-DIMENSIONAL STRUCTURES 

Jian-Ming Jin , John L. Volakis, and Jeffery D. Collins 
Radiation Laboratory 

Department of Electrical Engineering and Computer Science 
The University of Michigan 
Ann Arbor , Michigan 48109-2122 

ABSTRACT 

This paper presents a review of a hybrid finite element - boundary integral formu 
lation for scattering and radiation by two- and three-dimensional composite structures. 
In contrast to other hybrid techniques involving the finite element method, the proposed 
one is in principle exact and can be implemented using a low O(N) storage. This is of 
particular importance for large scale applications and is a characteristic of the boundary 
chosen to terminate the finite element mesh, usually as close to the structure as possible. 
A certain class of these boundaries lead to convolutional boundary integrals which can 
be evaluated via the fast Fourier transform (FFT) without a need to generate a matrix, 
thus, retaining the O(N) storage requirement. The paper begins with a general descrip- 
tion of the method. A number of two- and three-dimensional applications are then given, 
including numerical computations which demonstrate the method’s accuracy, efficiency 
and capability. 
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APPLICATIONS OF THE CONJUGATE GRADIENT FFT METHOD IN 
SCATTERING AND RADIATION INCLUDING SIMULATIONS 
WITH IMPEDANCE BOUNDARY CONDITIONS 

Abstract 


The theoretical and computational aspects related to the application of the 
Conjugate Gradient FFT (CGFFT) method in computational electromagnetics 
are examined. The advantages of applying the CGFFT method to a class of large 
scale scattering and radiation problems are outlined. The main advantages of 
the method stem from its iterative nature which eliminates a need to form the 
system matrix (thus reducing the computer memory allocation requirements) and 
guarantees convergence to the true solution in a finite number of steps. Moreover, 
since the CGFFT algorithm is highly vectorizable, it can be efficiently implemented 
on supercomputers and multiprocessor machines. 

Results are presented for various radiators and scatterers including thin cylin 
drical dipole antennas, thin conductive and resistive strips and plates, as well as 
dielectric cylinders. 

Solutions of integral equations derived on the basis of generalized impedance 
boundary conditions (GIBC) are also examined. These boundary conditions can 
be used to replace the profile of a material coating by an impedance sheet or insert, 
thus, eliminating the need to introduce unknown polarization currents within the 
volume of the layer. Moreover, by applying these surface boundary conditions, 
the difficulties associated with the calculation of the Green’s function are avoided. 
Impedance boundary conditions of up to the third order are employed and shown to 
be compatible with the basic CGFFT formulation, allowing an efficient simulation 
of large coated structures and filled cavity-backed apertures by further reducing 
the memory demand. For the purpose of validation of these simulations, a general 
full-wave analysis of two- and three-dimensional rectangular grooves and cavities 
is presented which will also serve as reference for future work. 
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Abstract 

Edge-based finite elements are free from the drawbacks associated with node- 
based vectorial finite elements and are. therefore, ideal for solving three- 
dimensional scattering problems. The finite element discretization using edge 
elements is checked by solving for the resonant frequencies of a closed in- 
homogeneously filled metallic cavity. Significant improvements in accuracy 
are observed when compared to the classical node-based approach with no 
penalty in terms of computational time and with the expected absence of 
‘spurious’ modes. A performance comparison between the edge-based tetra- 
hedra and rectangular brick elements is carried out and tetrahedral elements 
are found to be more accurate than rectangular bricks for a given storage 
intensity. A detailed formulation for the scattering problem with various 
approaches for terminating the finite element mesh is also presented. 
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